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INTRODUCTION 


The  linear  fracture  mechanics  problem  of  a  pressurized  cylinder  with  many 
radial  cracks  has  been  extensively  studied  in  recent  years  by  a  number  of 
investigators.  In  earlier  works,  cracks  were  assumed  to  be  uniformly 
distributed  and  of  equal  depths  so  the  problem  remained  axisymmetric.  In  a 
recent  report  (ref  1),  radial  cracks  of  unequal  depths  were  considered.  It 
was  found  that  the  variation  in  stress  intensity  factor  at  a  crack  tip  due  to 
crack  depth  change  of  the  same  crack  or  another  crack  is  approximately  linear. 
A  simple  numerical  method  was  used  in  Reference  1  to  accurately  estimate 
stress  intensity  factors  at  radial  cracks  of  unequal  depths  in  a  pressurized, 
non-autof rettaged  cylinder. 

For  a  partially  autof rettaged  cylinder,  the  functional  stress  intensity 
approach  was  developed  (ref  2)  for  equal  depth  radial  cracks.  The  same 
approach  can  be  extended  to  radial  cracks  of  unequal  depths.  The  functional 
stress  intensity  factors  for  uneven  cracks  must  be  first  obtained  by  the 
numerical  method  developed  in  Reference  1.  A  combination  of  methods  described 
in  References  1  and  2  can  be  used  to  estimate  stress  intensity  factors  at 
uneven  radial  cracks  in  a  pressurized,  partially  autof rettaged  cylinder. 
Numerical  results  are  given  for  two  diametrically  opposite  cracks  with  crack 
length  ratio  varying  from  0.7  to  1.4  for  a  thick-walled  cylinder  of  wall  ratio 
W  =  2.  The  method,  however,  is  general  for  a  multi-crack  configuration. 

Ipu,  S.  L.,  "Stress  Intensity  Factors  For  a  Circular  Ring  With  Uniform  Array 
of  Radial  Cracks  of  Unequal  Depth,"  ARDC  Technical  Report  No.  ARLCB-TR-84021 , 
Benet  Weapons  Laboratory,  Watervliet,  NY,  June  1984. 

2Pu,  S.  L.,  "A  Functional  Stress  Intensity  Approach  to  Multiply-Cracked, 
Partially  Autof rettaged  Cylinders,"  Transactions  of  the  28th  Conference  of 
Army  Mathematicians,  ARO  Report  83—1,  1983,  pp.  263—283. 
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RESIDUAL  STRESSES  IN  AN  AUTOFRETTAGED  CYLINDER 


It  is  now  a  routine  procedure  to  autofrettage  a  cannon  tube  to  introduce 
residual  compressive  stress  in  the  circumferential  direction  near  the  bore  to 
slow  the  growth  of  radial  cracks  which  occur  after  repeated  firing.  The 
residual  stress  distribution  is  difficult  to  measure  and  can  only  be  estimated 
by  theoretical  prediction  based  on  various  assumptions  of  material  models. 

From  a  recent  paper  (ref  3) ,  it  is  assumed  that  the  residual  stresses  can  be 
closely  approximated  by  the  following  expressions: 

1  <  r  <  1+et  (1) 

1+et  <  r  <  W  (2) 

1  <  r  <  1+et  (3) 

1+et  <  r  <  W  (4) 

where  (r,0)  are  polar  coordinates  with  origin  at  the  center  of  the  cylinder. 
The  inner  radius  is  taken  as  the  length  unit.  The  outer  radius  is  then  equal 
to  the  wall  ratio  W.  The  wall  thickness  t  equals  W-l.  The  degree  of 
autofrettage  e  is  defined  by  e  =  (rp-l)/t  where  rp  is  the  radius  of  elastic- 
plastic  interface  during  pressurization  in  an  autofrettage  procedure.  aQ  is 
the  yield  stress  of  the  material  and  or>  ag  are  stress  components  in  r  and  in¬ 
directions,  respectively.  The  superposition  coefficients  Ar,...,Er,  and 
Aq,...,Eq,  vary  with  W,  e,  and  the  material  model  used.  For  a  highly 
idealized  material  which  obeys  von  Mises'  criterion  and  is  incompressible  and 
elastic-perf ectly  plastic,  and  for  the  cylinder  which  is  under  a  state  of 

■^Pu,  S.  L.  and  Chen,  P.  C.  T. ,  "Stress  Intensity  Factors  For  Radial  Cracks  in 
a  Pre-stressed,  Thick-Walled  Cylinder  of  Strain-Hardening  Materials,"  J .  of 
Pressure  Vessel  Technology,  Vol.  105,  1983,  pp.  117-123. 


°r(r)/°o  = 


ae(r)/°o  = 


Ar  +  Br/r2  +  Cr!nr 
Dr  +  Er/r2 
Aq  +  Bg/r  2  +  C  g£nr 
D  0  +  Eg/r  2 


2 


plane  strain,  we  conclude  from  Reference  4  that: 


1 


2 


Aq  =  --  f 2  -  Pi(e01  ,  Bfl  =  —  Pi(e)  ,  Ce  =  -- 


(5) 


(6) 


where 


(1+et)  2 


+  2 £n(  1+et)  ] 


(7) 


For  a  less  restrictive  material  model  which  behaves  more  closely  to  high 
strength  steels  of  pressure  vessels,  residual  stresses  are  usually  obtained  by 
numerical  methods  and  are  given  by  numerical  values  at  discrete  points.  The 
superposition  coefficients  are  a  set  of  constants  found  by  the  least  squares 
method.  In  a  specific  example  (ref  3),  the  material  is  assumed  to  have  v  = 
0.3,  Et/E  =  0.1  where  Et  is  the  tangent-modulus  of  a  strain-hardening 
material,  and  the  material  follows  the  incremental  plastic  stress-strain 
relation  of  Prandtl-Reuss  and  the  yield  criterion  of  von  Mises  in  the  plastic 
region.  The  numerical  residual  stresses  obtained  by  Chen  (ref  5)  using  the 
finite  difference  method  for  a  cylinder  of  W  =  2  can  be  closely  approximated 
by  Aq  =  0.170,  Bq= 1.030,  and  Cq  =  0.890  for  e  =  1  and  by  a  different  set  of 
superposition  coefficients  for  a  different  value  of  e.  The  corresponding 
coefficients  from  Eqs.  (5)  and  (7)  for  idealized  materials  for  W  =  2,  e  =  1 

3Pu,  S.  L.  and  Chen,  P.  C.  T. ,  "Stress  Intensity  Factors  For  Radial  Cracks  in 
a  Pre-stressed,  Thick-Walled  Cylinder  of  Strain-Hardening  Materials,"  J .  of 
Pressure  Vessel  Technology,  Vol.  105,  1983,  pp.  117-123. 

^Hill,  R. ,  The  Mathematical  Theory  of  Plasticity,  Oxford  at  the  Clarendon 
Press,  1950. 

3Chen,  P.  C.  T. ,  "Numerical  Prediction  of  Residual  Stresses  in  an 
Autof rettaged  Tube  of  Compressible  Material,"  Proceedings  of  the  1981  Army 
Numerical  Analysis  and  Computer  Conference,  pp.  351-362. 
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are  Aq  =  0.0875,  Bq  =  -1.067,  and  Cq  =  1.155  which  are  quite  different  from 
the  coefficients  for  the  strain-hardening  material  given  previously. 

The  underlying  concept  here  is  that  the  residual  stresses  in  a  partially 
autof rettaged  cylinder,  without  reverse  yielding  during  elastic  unloading,  can 
be  closely  approximated  by  expressions  of  Eqs.  (1)  through  (4).  Once  the 
superposition  coefficients  are  found,  then  the  effect  of  such  residual 
stresses  on  stress  intensity  factors  at  radial  cracks  in  the  cylinder  can  be 
estimated . 


FUNCTIONAL  STRESS  INTENSITIES 

According  to  the  weight  function  method  (refs  6,7),  the  mode  I  stress 
intensity  factor  due  to  a  symmetrical  load  system  2  can  be  found  if  the  mode  I 
stress  intensity  factor  K*  and  displacement  field  u*  associated  with  the 
symmetric  load  system  1  for  the  same  cracked  geometry  are  known.  For  radially 
cracked  rings,  this  method  gives: 

H  ,a  3v* 

K  =  —  /  pc(x)  - —  dx  (8) 

K*  o  3a 

where  a  is  the  crack  length,  x  is  a  distance  measured  along  the  crack  from  the 
base  toward  the  tip,  v*  is  the  normal  component  of  crack  face  displacement,  H 
=  E  for  plane  stress  and  H  =  E/(l-v2)  for  plane  strain,  and  pc(x)  is  the 
crack  face  pressure  which  can  be  found  from  the  hoop  stress  in  an  uncracked 
ring  subjected  to  the  symmetric  load  of  interest. 

Taking  the  right  side  of  Eq.  (3)  as  pc(x)  into  Eq.  (8)  with  r  =  1+x,  we 


^Bueckner,  H.  F.,  "A  Novel  Principle  for  the  Computation  of  Stress  Intensity 
Factors,"  Z.  Agnew.  Math.  Mech.,  Vol.  50,  1970,  pp.  529-546. 

^Rice,  J.  R. ,  "Some  Remarks  on  Elastic  Crack-Tip  Stress  Fields,"  Int.  Journal 
of  Solids  and  Structures,  Vol.  8,  1972,  pp.  751-758. 
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have : 


----  =  AqK(1)  +  B  eK(r“ 2)  +  CqK( tor)  (9) 

°o 

where  K(e)  is  the  mode  I  stress  intensity  factor  at  a  radial  crack  emanating 
from  the  bore  with  crack  length  c  <  et  in  a  stress  field  given  by  the  residual 
stress,  Eq.  (3),  and  K(l) ,  K(r“2),  and  K( £nr)  are  integrals  given  by  Eq.  (8) 
with  pc(x)  =  1,  (1+x)-2  and  £n(l+x),  respectively.  These  integrals  are  stress 
intensity  factors  associated  with  the  specific  crack  face  loadings.  They  are 
basic  quantities  and  are  named  functional  stress  intensities. 

One  way  to  find  functional  stress  intensity  factors  is  to  assume  v*  as  a 
function  of  a  and  x  and  perform  the  integration  of  Eq.  (8).  An  approximation 
by  assuming  v*  as  a  conic  section  is  given  by  Grandt  (ref  8).  A  way  to 


deviate  from  this  approach  is  given  in  Reference  2  by  using  finite  element 
computations  of  stress  intensity  factors  of  the  cracked  cylinder  subjected  to 
the  following  three  loadings.  The  hook  stress  in  an  uncracked  cylinder  under 
internal  pressure  p^  is: 


°0  1  W2 

-  =  -2—  (1  +  -2) 
Pi  W2-l  r2 


(10) 


The  stress  intensity  factor  produced  by  this  stress  acting  on  the  crack  face 
is : 


K(Pi) 

~?i 


W' 


K(  1)  +  --—  K(r“2) 


W2-l 


W2-l 


(11) 


2Pu,  S.  L.,  "A  Functional  Stress  Intensity  Approach  to  Multiply-Cracked, 
Partially  Autof rettaged  Cylinders,”  Transactions  of  the  28th  Conference  of 
Army  Mathematicians,  ARO  Report  83—1,  1983,  pp.  263—283. 

8Grandt,  A.  F.,  "Stress  Intensity  Factors  for  Through-Cracked  Fastener  Holes, 
Int.  Journal  of  Fracture,  Vol.  11,  1975,  pp.  283-294. 
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Similarly,  the  stress  intensity  factor  due  to  a  uniform  tension  p0  on  the 
outer  cylindrical  surface  is: 


K(Po) 

Po 


W' 


U' 


K(l)  +  -3—  K(r 
W2-l  W2-l 


-  2\ 


(12) 


The  left-hand  sides  of  Eqs.  (11)  and  (12)  are  calculated  from  two  separate 
finite  element  computations.  The  functional  stress  intensity  factors  K(l)  and 
K(r-2)  are  solved  from  Eqs.  (11)  and  (12).  The  third  finite  element 
computation  which  supplies  values  to  the  left-hand  side  of  Eq.  (9)  is  obtained 
by  applying  an  axi symmetrical  thermal  loading  to  the  cylinder. 


T(r)  =  Tq  -  (TQ-Tw)*nJ)  1  <  r  <  W 


(13) 


where  TQ  is  the  temperature  at  the  bore  and  Tw  is  the  temperature  at  the  outer 
cylindrical  surface.  Temperature  gradient  is  related  to  material  constants 
by: 


To  -  TW  -  ~ 


4ao  (l-v)JtnW 
/3  Ea 


(14) 


where  E  and  a  are  Young* s  modulus  and  the  coefficient  of  linear  thermal 
expansion,  respectively-  The  thermal  stresses  produced  by  the  loading  given 
by  Eq.  (13)  are  identical  to  the  fully  autof rettaged  residual  stresses  in  Eqs. 
(1)  through  (4)  with  superposition  coefficients  given  by  Eqs.  (5)  and  (6);  see 
Reference  9.  The  only  unknown  K( £nr)  in  Eq.  (9)  can  be  easily  solved. 


^Hussain,  M.  A.,  Pu,  S.  L. ,  Vasilakis,  J.  D. ,  and  O’Hara,  P.,  “Simulation  of 
Partial  Autof rettage  by  Thermal  Loads,'*  Journal  of  Pressure  Vessel 
Technology,  Vol.  102,  No.  3,  1980,  pp.  314-318. 
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RADIAL  CRACKS  OF  UNEQUAL  DEPTHS 

For  non-autof rettaged  cylinders,  an  approximate  method  is  described  in 
Reference  1  to  estimate  stress  intensity  factors  at  radial  cracks  of  uneven 
depths.  Finite  element  code  APES  is  used  to  obtain  changes  in  stress 
intensity  factors  at  different  crack  tips  due  to  a  systematic  change  in  crack 
length  of  a  radial  crack.  The  total  change  in  stress  intensity  factor  at  a 
particular  crack  due  to  simultaneous  changes  in  crack  lengths  of  all  radial 
cracks  is  approximately  the  sum  of  changes  in  stress  intensity  factor  at  that 
crack  due  to  crack  length  changes,  one  at  a  time. 

For  partially  autof rettaged  cylinders  with  radial  cracks  of  equal  depth, 
we  use  functional  stress  intensities  defined  previously  as  basic  quantitites. 
The  final  stress  intensity  factor  is  given  by  a  linear  superposition  of 
functional  stress  intensities  such  as  Eq.  (9).  The  superposition  coefficients 
are  known  from  the  residual  stress  solution  and  the  Lame "  solution.  For 
radial  cracks  of  unequal  depths,  the  approximate  method  developed  in  Reference 
1  is  applied  to  each  functional  stress  intensity  factor.  Once  changes  in  all 
functional  stress  intensity  factors  are  known  due  to  simultaneous  changes  in 
crack  lengths  of  radial  cracks,  the  superposition  formula  can  be  used  to 
estimate  the  final  stress  intensity  factor  at  each  crack  tip  in  a  partially 
autof rettaged  and  pressurized  cylinder  with  uniform  array  of  radial  cracks  of 
unequal  depths. 

Ipu,  S.  L. ,  "Stress  Intensity  Factors  For  a  Circular  Ring  With  Uniform  Array 
of  Radial  Cracks  of  Unequal  Depth,”  ARDC  Technical  Report  No.  ARLCB-TR-84021 , 
Benet  Weapons  Laboratory,  Watervliet,  NY,  June  1984. 
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NUMERICAL  RESULTS  AND  DISCUSSIONS 

Using  two  diametrically  opposite  cracks  as  an  example,  the  stress 
intensity  factor  is  plotted  as  a  function  of  crack  depth  in  Figure  1  for  an 
equal  crack  configuration  of  a  non— autof rettaged  cylinder  of  W  =  2  subjected 
to  a  uniform  tension  p  on  the  outer  cylindrical  surface.  In  the  same  figure, 
results  of  n  -  1,  3,  4,  6  are  also  shown  to- indicate  that  stress  intensity 
factor  for  n  =  2  is  the  highest  at  a  given  crack  depth  c.  The  subscript  e  is 
used  to  emphasize  that  the  quantity  is  of  equal  crack  depth  configuration.  In 
Figure  2,  the  change  in  stress  intensity  ratio  is  plotted  as  a  function  of 
change  in  crack  depth  when  one  of  the  cracks  is  growing  from  the  equal  depth 
to  an  unequal  depth  configuration.  To  distinguish  one  crack  from  other  cracks 
in  an  unequal  depth  configuration,  an  arbitrary  numbering  system  is  used. 
Quantities  with  a  subscript  1,  such  as  and  c^ ,  are  the  ones  associated  with 
crack  1.  Dimensionless  stress  intensity  factor  and  crack  depth  are  defined  by 
Ni  =  Kj_/Ke  and  =  c^/c.  The  changes  in  N  and  p  are  designated  by  AN^  =  (K^- 
Ke)/Ke  anci  APi  =  (ci-c)/c. 

For  partially  autof rettaged  cylinders,  functional  stress  intensity 
factors  are  studied  as  basic  quantities.  Figure  3  shows  functional  stress 
intensity  factors  as  functions  of  crack  depth  when  two  radial  cracks  are  of 
equal  depths.  Results  in  Figure  3  are  used  to  normalize  functional  stress 
intensities  when  radial  cracks  are  of  unequal  dept*hs.  Figures  4  through  6, 
similar  to  Figure  2  of  non-autof rettaged  cylinders,  are  graphs  of  ANj_  versus 
Api  for  the  three  basic  functional  stress  intensities  corresponding  to  normal 
crack  face  loadings  pc  =  1,  1/r2,  and  £nr,  respectively. 
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Using  the  approximate  method  described  in  Reference  1  and  the  results  in 
Figures  4  through  6,  functional  stress  intensity  factors  can  be  estimated  for 
unequal  depth  radial  cracks*  The  approximate  results  agree  well  with  results 
directly  calculated  from  finite  element  computations*  For  instance,  for  n  — 

2 ,  c2  =  c  =  0.2  and  ci  varies  from  0.14  to  0.28.  A  comparison  of  approximate 
and  finite  element  results  of  functional  stress  intensities  is  given  in  Table 


I. 


TABLE  I.  APPROXIMATE  VS.  FINITE  ELEMENT  RESULTS  OF  K^p^/KeCp,.) 


FOR  W  =  2,  n  =  2,  c  =  0.2,  =  0.7  and  =  1.4 


i 

1 

1 

1 

i 

C2  =  0.2,  pi  =  0.7  | 

C2  -  0*2, 

Pi 

=  1.4 

1 

pc 

Approximate 

i  r 

|  F.E.  I 

Approximate 

1 

1 

F.E. 

T 

l  1 

| 

0.8113 

i  r 

1  0.8130  | 

1.2515 

1 

1 

| 

1.2519 

i 

r"2  I 

0.8734 

1  1 
|  0.8622  | 

1.1688 

1 

1 

1 

1.1653 

inv  | 

0.5082 

1  1 

I  0.5698  | 

1.6557 

1 

1 

1.6797 

1  | 
i 

0.9629 

I  0.9710  | 

1  | 

1.0495 

1 

1 

1.0519 

2 

1 

r-2  I 

1 

0.9643 

1  1 
|  0.9706  I 

1 

1.0477 

1 

1 

1 

1.0495 

1 

Inr  | 

1 

0.9574 

1  1 

|  0.9473  | 

1  1 

1.0568 

1 

1 

1 

1.0617 

L 


Knowing  Ki(pc)/Ke(pc)  and*Ke(pc),  Figure  3,  then  Eq.  (9)  can  be  used  to 
compute  stress  intensity  factors  at  radial  cracks  of  unequal  depths  corres¬ 
ponding  to  a  residual  stress  distribution  of  Eq.  (3).  The  superposition 
coefficients  are  known  from  the  residual  stress  distribution.  For  the 


1-Pu,  S.  L. ,  "Stress  Intensity  Factors  For  a  Circular  Ring  With  Uniform  Array 
of  Radial  Cracks  of  Unequal  Depth,"  ARDC  Technical  Report  No.  ARLCB-TR-84021 , 
Benet  Weapons  Laboratory,  Watervliet,  NY,  June  1984. 
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idealized  material  model,  they  are  given  by  Eq.  (5).  Let  e  =  1.0  and  0.6; 
then  numerical  values  of  K^( e)/ ( 0o/Rq)  are  obtained  by  the  approximate  method 
and  the  finite  element  method.  The  results  are  compared  in  Table  II.  The 
negative  sign  indicates  that  the  effect  of  autofrettage  residual  stress  on 
inner  cracks  is  to  resist  crack  opening  due  to  internal  pressure  or  other 
loadings.  The  resistance  varies  and  never  .exceeds  the  crack  opening  force. 
The  values  in  Table  II  are  the  maximum  resistance.  The  stress  intensity 
factor  at  a  crack  remains  zero  unless  the  crack  opening  force  is  larger  than 
the  maximum  resistance  given  by  the  residual  stress  distribution. 


TABLE  II.  APPROXIMATE  VS.  FINITE  ELEMENT  RESULTS  OF  Kt( e)/( OqRi1/2) 
FOR  W  =  2,  a  =  2,  c  =  0.2,  px  =  0.7  and  Pi  =  1.4 


C2  -  0.2,  P]_  =  0.7 


Approximate 


F.E< 


C2  =  0.2,  =  1.4 


Approximate 


F.E. 


1 

1 

1 

1 

1.0  | 

1 

-0.6201 

1  1 

I  -0.6001  | 

1  | 

-0.6938 

1  1 
|  -0.6879  | 

1 

1 

1  1 

1 

1 

0.6  | 

i 

1 

-0.5236 

I  -0.5040  | 

1  1 

-0.5545 

1 

I  -0.5486  | 

1  1 

- 1 

1 

1 

—  i . 

1.0  | 

1 

-0.6294 

1  | j 

|  -0.6320  | 

1  | 

-0.6814 

1  — — | 

|  -0.6822  | 

1 

1 

2  I 

1 

0.6  | 

1 

-0.5187 

I  -0.5205  | 

1  1 

-0.5610 

1 

|  -0.5615  | 

1  1 

For  an  inner  cracked,  partially  autof rettaged  cylinder  subject  to  an 
internal  pressure  p,  the  resultant  stress  intensity  factor  is  the  sum  of  Eqs. 
(9)  and  (12).  If  p  Is  oQ/ f  where  f  is  a  constant,  the  resultant  stress 
intensity  factor  at  crack  i  is  given  in  terms  of  functional  stress  intensity 
factors  by: 
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Ki(p,e) 


Ki(l) 


- -  (A0(e.)  + - t - ) - +  (B  e(  e)  + - - - ) 

a0/ST  f(W2-D  /r^  f  (W2— 1) 


KiCr-2) 

- +  C  q(  e) 


KiC  £ar) 

~’^T 


(15) 


Normalizing  K.£(p,e)  by  KgCpj^O),  numerical  results  for  i  =  1  or  2,  e  =  0.5  or 
1,  and  f  -  1.5  or  3  are  shown  in  Figures  7  through  9  for  c  -  0.1,  0.2,  and 
0.3,  respectively.  In  these  figures,  the  superposition  coefficients  Aq,  Bq, 
and  Cq  are  given  by  Eq.  (5).  For  a  strain-hardening  material  with  tangent 
modulus  Et  =  0.1E,  the  residual  stress  distributions  obtained  by  Chen  (ref  5) 
were  closely  approximated  by  various  sets  of  superposition  coefficients  in 
Reference  3.  Using  those,  the  final  stress  intensity  factors  ratios  are  shown 
in  solid  lines  in  Figure  10  in  comparison  with  broken  lines  for  idealized 
materials . 

In  Figures  7  through  9,  all  lines  with  the  same  i  are  nearly  parallel, 
which  indicates  that  the  effect  of  autofrettage  residual  stress  is  simply  a 
proportional  reduction  of  stress  intensity  factor.  The  same  trend  is  true  in 
Figure  10  for  strain-hardening  material.  The  resistance  to  crack  opening  due 
to  autofrettage  is  slightly  reduced  in  strain-hardening  material  than  in  the 
idealized  material  (ref  3) .  This  reduction  is  in  proportion  to  the  change  in 
crack  depths  so  the  solid  lines  remain  parallel  to  the  broken  lines  in  Figure 
10. 


2 

Pu,  S.  L.  and  Chen,  P.  C.  T.,  "Stress  Intensity  Factors  For  Radial  Cracks  in 
a  Pre-s tressed ,  Thick— Walled  Cylinder  of  Strain-Hardening  Materials,"  J.  of 
Pressure  Vessel  Technology,  Vol.  105,  1983,  pp.  117-123. 

^ Chen,  p"!  C .  T. ,  "Numerical  Prediction  of  Residual  Stresses  in  an 
Autofrettaged  Tube  of  Compressible  Material,"  Proceedings  of  the  1981  Army 
Numerical  Analysis  and  Computer  Conference,  pp.  351-362. 
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Figure  1.  Stress  intensity  factors  as  a  function  of  crack  depth  for  a 

pressurized,  non-autofrettaged  cylinder  of  W  =  2  with  n  cracks 
of  equal  depths.  Dots  are  finite  element  results. 
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ANi<p)  = 


Kj(p)  -  Ke(p) 
Ke(P) 


Figure  2.  Changes  in  stress  intensity  ratio  as  a  function  of 
relative  change  in  one  of  the  crack  depths  for 
pressurized  cylinder  with  two  diametrically  opposite 
cracks.  Dots  are  finite  element  results. 
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Figure  3.  Functional  stress  intensity  factors  as  a  function  of 
crack  depth  for  two  equal  cracks  under  three  types  of 
basic  crack  face  loadings.  Dots  are  finite  element 
results . 
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ANjIPc  > 


Figure  4,  Changes  in  stress  intensity  ratio  as  a  function  of 

relative  change  in  one  of  the  crack  depths  for  a  two- 
crack  configuration  due  to  a  constant  crack  face 
"*  loading.  Dots  are  finite  element  results. 
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Figure  5.  Changes  in  stress  intensity  ratio  as  a  function  .of 

relative  change  in  one  of  the  crack  depths  for  a  two- 
crack  configuration  due  to  a  crack  face  loading  of  pr-2. 
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Figure  6.  Changes  in  stress  intensity  ratio  as  a  function  of 

relative  change  in  one  of  the  crack  depths  for  a  two- 
crack  configuration  due  to  a  crack  face  loading  of 
pln(r) . 
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Kilp-oyf.C) 


F-igure  7.  Normalized  stress  intensity  factors  as  a  function  of 
relative  change  in  one  of  the  crack  depths,  with  the 
other  crack  fixed  at  a  depth  c  =  0.1,  in  a  diametrically 
opposite  cracked  cylinder  which  is  autofrettaged  to  a 
degree  of  e  and  is  subjected  to  an  internal  pressure 

P  =  Oq/ f . 
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Kj(p=Ob/f,£) 

Ke(p,£=0) 


Figure  8.  Normalized  stress  intensity  factors  as  a  function  of 
relative  change  in  one  of  the  crack  depths,  with  the 
other  crack  fixed  at  a  depth  c  *  0.2,  in  a  diametrically 
opposite  cracked  cylinder  which  is  autofrettaged  to  a 
degree  of  e  and  is  subjected  to  an  internal  pressure 
P  ■  crQ/f. 
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Ki(p=cr0/f,£) 


Figure  9.  Normalized  stress  intensity  factors  as  a  function  of 
relative  change  in  one  of  the  crack  depths,  with  the 
other  crack  fixed  at  a  depth  c  =  0.3,  in  a  diametrically 
opposite  cracked  cylinder  which  is  autofrettaged  to  a 
degree  of  £  and  is  subjected  to  an  internal  pressure 
P  =  °0/f. 
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Figure  10. 
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Comparison  of  normalized  stress  intensity  factors  for  strain 
hardening  materials  having  Et  -  0.1E  with  corresponding 
values  for  elastic— perfectly  plastic  materials. 
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